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Abstract 

Considering the physical 3-space t = constant of the space-time 
metrics as spheroidal and pseudo spheroidal, cosmological models which 
are generalizations of Robertson- Walker models are obtained. Specific 
forms of these general models as solutions of Einstein's field equations 
are also discussed in the radiation and the matter dominated era of 
the universe. 



1 Introduction 

In general relativity, a gravitational situation is equivalent to imposing a 
Riemannian structure on the space-time fourfold. There exists a correspon- 
dence between matter-energy and the geometry of the physical 3-space t 
= constant. The space-times associated with the Robert son- Walker metrics, 
representing standard models of the universe, contain matter-energy to curve 
up the physical 3-space t=constant into 3-spheres and 3-pseudo spheres [I]. 
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The absence of the dissipative terms in the energy-momentum tensor of the 
standard models, due to the assumption that the matter in the universe is ho- 
mogeneous and isotropic, suggests that these models may not appropriately 
describe the early hot dense state of the big bang cosmology j2] . Cosmological 
models described by the space-times which have hypersurfaces t=constant as 
3-spheroids and 3-pseudo spheroids are obtained in order to incorporate the 
dissipative terms in the energy-momentum tensor. Robert son- Walker models 
are special cases of these general models. 



2 Derivation of the metrics 

A 3-spheroid, immersed in a 4-dimensional Euclidean flat space given by the 
metric 



{dS E f = (dx 1 ) 2 + (dx 2 ) 2 + (dx 6 ) 2 + (dx 4 ) 2 



(2.1) 



in Cartesian coordinates ( ), will have the Cartesian equation 

,.4\2 



(x 1 ) 2 + (x 2 ) 2 + (x^ 2 



+ 



[X 



(2.2) 



R 2 S 2 

where R and S are constants. We use coordinates intrinsic to the surface of 
the 3-spheroid 

,i 



x 



x 



X 

and x 4 



.Rsinxcos^, 
R sin x sin 6 cos < 
R sin x sin 9 sin ^ 
S cos x- 



(2.3) 



The spatial line-element on the surface of the 3-spheroid is obtained from the 
metric (|2.1j) of the 4-dimensional flat space in which it is immersed as 



(dSi 



R 2 



S 2 

cos 2 x + s ^ n2 ^ ( ^ 2 s ^ n2 ^(^ 2 + s ^ n2 ° ' ? 1 



(2.4) 
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We now write sin^ = r(0 < r < 1) and obtain (cLSe) 2 in r, 9, <j) coordinates 
as 



) 2 = R 2 



1 + 



ff 2 r 2 



dr 2 + r 2 (d9 2 + sin 



(2.5) 



The space-time metric, which has its associated physical 3-space t=constant 
3-spheroid, is therefore written as 



ds 2 = -dt 2 + R 2 



1 + 



^2 r 2 



dr 2 + r 2 (d8 2 + sin 



(2.6) 



R 2 1 - r 2 

where i? and S 1 are now functions of time t. Relativistic units are being used 
with c — G—1. 

Next we consider a 3-pseudo spheroid, with the Cartesian equation 
(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 (x 4 ) 2 



-1 



(2.7) 



R 2 S 2 

where R and S are constants, immersed in a 4-dimensional pseudo Euclidean 
flat space with metric 



(dS p ) 2 = (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 - (dx 4 ) 2 
We use coordinates intrinsic to the surface of the 3-pseudo spheroid 



(2.8) 



.x 



X 



X' 

and x 4 



R sin hxcos9, 
R sin hx sin 9 cos < 
R sin hx sin 9 sin c 
S cos hx- 



(2.9) 



The spatial line-element on the surface of the 3-pseudo spheroid is obtained 
from the metric ()2.8|) of the 4-dimensional pseudo Euclidean flat space in 
which it is immersed as 



(dSp) 



R 2 



S 2 

cos h 2 x — iTTr sin h 2 x f dx 2 + sin h 2 x(d9 2 + sin 2 
R l 



(2.10) 
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We now write sin hx = t and obtain (dS p ) 2 in r, 6, 4> coordinates as 

S 2 r 2 



(dS f 



) 2 = R 2 



dr 2 + r 2 (d9 2 + sin" 



(2.11) 



R 2 1 + r 2 

The space-time metric, which has its associated physical 3-space £=constant 
as a 3-pseudo spheroid, is therefore written as 

S 2 r 2 



ds 2 



-dt 2 + R 2 



dr 2 + r 2 (d9 2 + sin 2 



(2.12) 



R 2 1 + r 2 

where R and S are now functions of time t. The two line-elements (|2.6|) and 
(I2.12J) can however, be combined into a single line-element with the help of 
a parameter k that takes values k = ±1 as 

S 2 kr 2 



ds' 



-dt 2 + R 2 



dr 2 + r 2 (d9 2 + sin 2 



(2.13) 



^(l-Ar 2 ) 

It may be noted that if the parameter k is set equal to zero we get the 
line-element f!2.13|) as 

2 - -dt 2 + R 2 [dr 2 + r 2 (d9 2 + sin 2 6d<f) 2 )} (2.14) 



ds 2 



which has its associated physical 3-space t = constant as a 3-surface of zero 
curvature. The space-time metrics ()2.13|) in normal Gaussian coordinates 
r, 9, (ft, t are the natural generalizations of the Robertson- Walker space-times 
for which S(t) = R(t). The signature of the space-time line-element ()2.13j) 
is (+, +, +, -) as (1 + ^ tjz^) ) is always positive. 

3 The Field Equations 

The energy-momentum tensor of a relativistic imperfect fluid is given by [2] 
T ij = pg ij + {p + p^u 3 - r]h ik h jl w kl - x{h ik u j + h? k ij)q k - £h i3 u k ; k (3.1) 
where Wij is the shear tensor 



— 4- k - h 



(3.2) 



qi is the heat flow vector 

qi = T.i + Tui-jV? , (3.3) 

hij is the projection tensor on the hypersurface normal to u % 

hij = g { j + UiUj, (3.4) 

X is the coefficient of heat conduction, 77 is the coefficient of shear viscosity, 
£ is the coefficient of bulk viscosity, T is the temperature, p is the isotropic 
fluid pressure, p is the proper matter-energy density and u % is the fluid flow 
vector satisfying gijU l v? — — 1. A semi-colon stands for the covariant differ- 
entiation. The thermodynamic behaviour of a dissipative fluid is governed 
by the following set of thermodynamic laws [2] 

1. Particle number conservation: 

(nu% = (3.5) 

2. The second law of thermodynamics: 

KTda = pd (^j + d (3.6) 

where aK is the entropy per particle. K is the Boltzmann constant intro- 
duced here to make a dimensionless. An entropy currnet four-vector S % is 
defined by the equation 

S i = nKou 1 - xT~y . (3.7) 

The rate of entropy production per unit volume is given by S^. It may be 
noted that S t i is non-negative if x ^ 0- The models are also required to 
satisfy various physical conditions such as < p < | and r] > 0, £ > 0. 
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In the adopted comoving coordinate system we have the components of 
the fluid 4-velocity vector u 1 as u l = u 2 = u 3 = 0, u A = 1 and u^v? = 0. The 
volume expansion 9 = u l ;i for the models (j2.13|) is obtained as 

3R / S 2 kr 2 \~ l ( S R \ S 2 kr 2 



6 R + \ 1 + R 2 (l-kr 2 )) is RJ R 2 (1 - kr 2 ) ^ 

where an overhead dot denotes differentiation with respect to time t. The 
non-zero components of the shear tensor Wij for the models (j2.13j) are ob- 
tained as 

Wl ~ 3\ + R 2 {l-kr^) \S~ RJ R 2 (l-kri) ( } 

and 

, „ if S 2 kr 2 Y 1 ( S k\ S 2 kr 2 , . 

W2 = w * = -3{ 1+ i?jr^)) {s-r)r- 2 Jti^ A0) 

From equation (|3.5|) . the particle number density n for the models (j2.13|) is 
obtained as 

/ C2 u 2 \ -1/2 

" = ^ 3 ( 1 + ^(T^J ' 3,11) 

where A(r) is an arbitrary function of r. Since the models ()2.13|) reduce to 
R — W models when S = R we obtain A(r) = (1 — kr 2 )^ 1 ^ 2 and therefore, 

„ = fl -3 (1 -^)-V,( 1 + |_^_)- 1/2 (3 , 2) 

The non-vanishing components of the Einstein field equations 



- -Rg) = -8tt3] (3.13) 



for the line-elements fl2.13j) and the imperfect fluid distribution (|3.1j) are 
obtained as 

2R R 2 S 2 k ( S 2 kr 2 Y 1 . A/1 . 
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2R &_ S 2 kr 2 ( S 2 kr 2 \~ 1 [ S R R 2 RS\ 
R + R 2 + R 2 (1 - kr 2 ) V + R 2 (1 - kr 2 ) ) [S R R 2 + RS J 

/ S 2 kr 2 y 2 [ S 2 k S 2 kr 2 S R 2 \ 

+ V + R 2 (l-kr 2 )) jr 4 (l-£x 2 ) 2 + R 2 1 - kr 2 ^ S R ] } 

= -8n[p -£0 - rjwj], 

(3.15) 



3R 2 1 

+ 



R 2 R 2 r 2 



+ 1 + 



S 2 



kr 2 



R 2 (1 - kr 2 ) 



-i 



R 2 r 2 
2S 2 k 



S_R \ 2S 2 R kr- 



R\l - kr 2 ) 2 



S R R 2 R{1 

S 2 kr 2 



kr 2 ) 



1 + 



R 2 (1 - kr 2 ) 



and 



S 2 kr 2 



-i 



S 
S 



R\ 2S 2 
R 



kr 



R 2 (l-kr 2 )J \S RlR 2 (l-kr 2 ) 
From equations ()3.14j) and ()3.15j) we get 



(3.16) 



(3.17) 



S 2 kr 2 



„ . r „^ QR 3R 2 



-1 



S 2 k 



2S 2 kr 2 



+ 1 + 



S 2 kr 2 
R^{l-kr 2 ) 



R 4 (l - kr 2 ) 

2 



S R R 2 R S 
f R 2 (1 - kr 2 ) \S~R~R 2 + RS 



2S 2 k 



+ 



2S 2 kr 2 



S R 



R\l - kr 2 ) 2 R 2 (1 - kr 2 ) \S R 



(3.18) 
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and 

RJ ~ V + R 2 (l-kr 2 )) } R 2 r 2 (l-kr 2 ) + yS R 

S R R 2 R S 1 

+ 




S R R 2 RS \ R 2 r 2 

(3.19) 



From equations (|3.17|) and (|3.3j) we get 



1 S 2 kr 2 Y l ( S R\2S 2 kr o dT . . 

1 + f^(i^)J " i?J 1p (T3fo5) = 87rx a7- (3 ' 20) 

Integrating this equation we get 

.V /? \ ifi) , f / S 2 kr 2 \- l fS 2 



8 ^ T = [s-r) (f=ty log i 1 + l 1+ # (i^y ) U " 1 1 Kfl(t) 

(3.21) 

where -B(t) is an arbitrary function of time t. Since the models ()2.13|) reduce 
to R-W models when S = R we get B(t) = 8itxTr-w an d equation ()3.21|) is 
now expressed as 

**<T - 7W) = (| - |) ifO_ lo g { {l _ff_ g)} > ■ (3.22) 

4 Radiation Dominated Models 

If the energy-density of the universe is dominated by highly relativistic par- 
ticles then p = \p and £ = [2|. From equations ()3.16|) and ()3.18|) in this 



case we get after some simplification 



2„4 



k r 




2R S R 2 2RS 
~R + S + R 2 + ~RS 



+ 



R 2 



5R 5R 2 S 
~R + ~R 2 ~ + S 



2R 2 



3R 3R 2 
~R + ~R? 



S 2 



3R 3R 2 
~R + ~R 2 



+ kr z 




5R 5R 2 S S 2 
~R + ~R 2 ~ + S + S 2 



+ 



S 2 



3R 3R 2 
~R + ~R 2 



+ 



3A; 
R 2 



(4.1) 



This implies 

S 2 I 2/? 



R 2 



R? 



S R 2 2R S 
R + S + R 2+ ~RS | + S 2 



3R 3R 2 
~R + ~R? 



5R 5R 2 S S 2 
~R + l¥ + S + ~S 2 



(4.2) 



\~R + ~& + S + S 2 j + ~R^J ~R 2+ ~R r ~ ( ' 

and 

# 2 \3R 3R 2 \ 3k 

^{-R + ^) + R 2 = - (44) 

Equation (|4.4j) gives 

As expected ^{R 2 ) is positive or negative according as k — — 1 or +1. From 
equations (|4.2|) and ()4.3|) we get 

2R R 2 k S 2SR 2k n lA ^ 

lf + i? + 7? + s + TS + ^ = ' (4 ' 6) 

If we take § = , that is if we take S = aR(0 < a < 1) then equations ()4.5|) 
and ()4.6j) imply that a — 1, that is the models are R-W models. Making use 
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of equation (|4.5|) in (|4.6|) we obtain a differential equation for R 2 as 



+ !!<« 2 >> 2 > 



+ 



2k r d 2 



R 2 I eft 2 

Equation (|3.6J) in this case takes the form 



(R 2 ) - 



This implies 



and 



d(Ka) = ^d(- 



T = 3n 1/3 



nT 



dp. 



Kg = -^pn 



-4/3 



In view of equation ()3.12|) we get equation ()4.9|) 

S 2 kr 2 



T = 3 



E (1 - fcr 2 ) 1 + 



R 2 (1 - £x 2 ) 



1/6" 



(4.7) 



(4.8) 

(4.9) 
(4.10) 



(4.11) 



where S is given in terms of R from equation ()4.5j) and i? is given by the 
differential equation (|4.7J) . Making use of equation ()4.11|) in ()3.22|) we obtain 



iTrxlsR-^l-kr 2 )- 1 ^ ( 1 + 



5 R 



S 2 kr 2 
R 2 (1 - At 2 



-1/6 



3(i? r 



log 



i? 2 



(4.12) 



where iiw is the scale factor of the R-W models in radiation dominated era. 
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With the help of equations (|3.12|) and (|3.16|) we can write equation (j4.1U|) as 

a kr 2 x -i 



R 2 r 2 



Ka 



S_R \ 2S 2 R kr 2 



24tt 



3R 2 1 / N 
2S 2 k 



1 + 



S' 



R 2 (1 - kr 2 ] f 
- kr 2 x 



S R J R 2 R(l - kr 2 ) \ R\l - kr 2 ) 2 

R *(l _ kr 2 ) 2 /^ ( 1 + 



R 2 (1 - kr 2 ] f 

S 2 kr 2 \ 2/3 ' 
R^{1 -kr 2 ) ) 

(4.13) 

where S is given in terms of R from equation ()4.5|) and f? is given by the 
differential equation (|4.7j) . 

5 Matter Dominated Models 

If the energy density of the universe is dominated by the non-relativistic 
matter with negligible pressure and bulk viscosity then equation ()3.18|) gives 
after some simplification 



2S 2 
S 2 



2„4 



k r 





AR 2S R 2 2RS 
~R + T + ^ + ~R S 



+ kr 



-2 



ok m 2 

~R + 





2S 6R 2 2S 2 



GR 3_R^ 3kS 2 
+ + ~R 2+ R 4 




(5.1) 
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This implies 




AR 2S R 2 2R S 1 S 2 I 10 R 2S 
~R + T + ^R2 +^ sj 'R 2 j~R~ + T" 

6i? 2 2S 2 2RS) \§R 3R 2 



(5.2) 




R 2 S 2 R S R R? 



10R 2S 6R 2 2S 2 _ 2R S 1 J 6R 3R 2 

/W 2 

+ ■ 



R & R 4 



(5.3) 



and 



6# 3R 2 3kS 2 n ,„ 
^ + ^ + ^ = 0- (5-4) 



From equation (|5.4|) we get 

S 2 = -^(2RR + R 2 ) (5.5) 

This shows that (2RR + R 2 ) is positive or negative according as k is —1 or 
+1. From equations (|5.2j) and ()5.3|) we get 

4fl i? 2 k\ 2S 2RS 2k . . 

-R + R- 2 + R- 2 ) + ^ + JlS + S- 2=0 (5 ' 6) 

If we take f = ^ then in this case also equations ()5.5|) and ()5.6|) imply that 
the models are R — W models. If we substitute for S from equation ()5.5j) in 
(J5.6D we again get a differential equation of fourth order for R. 
Equation ()3.6|) in this case leads to 

T = constant (say M) (5.7) 

and 



(5.8) 



Equation (|3.22|) in this case gives 



8tv X (M - T R . 



w. 



S 
S 



R 



RJ {1 



s 2 \ 
Wf 



log 



{l-^(l-f)} 



(5.9) 



where Tr_vi/ is the constant temperature of the R-W models in the matter 
dominated era. Making use of equations ()3.12j) and (jH.lfij) in equation ()5.8j) 

we get 

-l 



Ka 



3R Z 



+ 



R 2 R 2 r 2 



+ 1 



S 2 kr 2 
~R?(l- kr 2 



R 2 r< 



+ 



S R\2S 2 R kr 2 
S~Rj ~WR(l-kr 2 ) 



+ 



2S 2 k 



1 + 



S 2 kr 2 



i? 4 (l - kr 2 ) 2 \ R 2 (-kr 2 

/ C2 u 2 \ 1/2 



(5.10) 



R 2 (1 - kr 2 ) 



In the above two equations ()5.9|) and ([5.10)1 . S is given in terms of R from 
equation ()5.5|) and R is given by the differential equation obtained from (|5.6|) 
after substituting for S from (j5.5|) . 



6 The Hubble and Deacceleration parame- 
ters 

A representative length I which represents the volume behaviour of the fluid 
is defined by the equation 

I* liU* 1„ I / 

7 sJ T = s s = r < 6J) 

For the models ()2.13|) I is obtained as 

' = c < r K 1 + ^(T^)) (62) 
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where C(r) is an arbitrary function of r. Since the models (|2.13|) reduce to 
R-W models when S = R we get C(r) = (1 — Ax 2 ) 1 / 6 and therefore, 



I = R{1- kr 2 ) 1/6 ( 1 + 



S 2 kr 2 
R?{1- kr 2 ) 

The Hubble parameter H is defined by the equation 



1/6 



H 



For the models (j2.13|) . H is given by 



I* 

7 



s 2 



kr 2 



R 2 (1 - kr 2 ) 



S R \ S 2 



kr 2 



S R R 2 {1- kr 2 ) ' 



The deacceleration parameter q is defined by the equation 

/ / \ / / \~ 

Q = ~ 



I 



For the models (j2.13|) q is obtained as 



, R IS 2 



kr 2 



R 3R 2 (1- kr 2 ) 



S 2 



kr 2 





R 2 \ 


1(1- 


R 2 J 



S R R 2 
S~R~~S 2 + R 2 



+ 2 



R 2 (1 - kr 2 

1 S 2 kr 2 
3^(1 - kr 2 ) 



-i 



S 2 



S_R 
S R 

kr 2 



1 + 

S 2 kr 
R 2 (1 -kr 2 ) 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



-2 



-1 



R 2 (1 - kr 2 ) 

2 h./yT.I ^ ' 



From equations (j6.1j) . ()6.4|1 and ()6.6|) we find that 



(6.7) 



(6.8) 



With the help of equations ()6.3j) . ()6.5|) and ()6.7|) the expression for matter- 
energy density given by equation (j3.16|) and the expression for the isotropic 
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pressure given by equation (|3.18|) can be expressed as 



R 2 \ I „ ( R\ 1 \ f (I 6 - R 6 (l-kr 2 ^ 



lirp = 3 — \\2H - - 1 + 1 + 



„2 



R 2 j \ \ Rj I 1 I 6 I I Z 6 i? 2 r 

(6.9) 
or 

8np = ^Hr-wY {2H{H R ^ - 1} + |l + _ j | ^ i| 

(6.10) 

and 

/ ^ / , TT 2 3i? 2 {/ 6 - R\l - kr 2 )} / 2.R 6 
8vr(p - ^ = 2(1 + g)tf 2 - — - ± ^ 1 + 



R? PR 2 r 2 V I 6 

6(#-|) 2 / 6 |\ 2i? 6 (l - kr 2 ) 1 „ /„ i?^ R 



1 -6 fi 



{/ 6 - i? 6 (l - £x 2 )} [ / 6 J \ i? 

(6.11) 

or 

M , 2 x2 {/ 6 - i? 6 (l - fcr 2 )} / 2i? 6 
8yr(p - £0) = 2(1 + q)H 2 - 3(H R ^ W ) 2 - I ^ ^ (^1 + — 

6(jy - H R - W ) 2 l G L 2R 6 (l-kr 2 



{/ 6 - /? 6 (1 - At 2 )} 



I }> - 6(ff - H R _ W )H R _ W . 

(6.12) 
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